A non inflationary model with scale invariant cosmological perturbations 
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We show that a contracting universe which bounces due to quantum cosmological effects and 
connects to the hot big-bang expansion phase, can produce an almost scale invariant spectrum of 
perturbations provided the perturbations are produced during an almost matter dominated era in 
the contraction phase. This is achieved using Bohmian solutions of the canonical Wheeler-de Witt 
equation, thus treating both the background and the perturbations in a fully quantum manner. We 
find a very slightly blue spectrum (n s — 1 > 0). Taking into account the spectral index constraint 
as well as the CMB normalization measure yields an equation of state that should be less than 
w<8x fO -4 , implying n s — f ~ O (lO -4 ), and that the characteristic size of the Universe at the 
bounce is Lq ~ 10 £ P1 , a region where one expects that the Wheeler-DeWitt equation should be 
valid without being spoiled by string or loop quantum gravity effects. 
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I. INTRODUCTION 

When the theory of cosmological perturbations Q] is 
applied to a background cosmological model and their 
initial spectra can be justified in physical terms, as it 
is the case of inflationary and nonsingular models where 
the horizon problem is solved, one obtains definite pre- 
dictions concerning the spectrum of primordial scalar 
and tensor perturbations in this background, establish- 
ing the initial conditions for structure formation and the 
angular power spectrum of the anisotropies of the cosmic 
microwave background radiation (CMBR). Since the re- 
lease of the first data obtained from observations of these 
anisotropies in 1992 pj, primordial cosmological models 
have been confronted with new observational facts 0, 
bringing them definitely out of the arena of speculations 
based almost uniquely on theoretical and aesthetical ar- 
guments. In fact, many interesting cosmological back- 
grounds have been falsified since then |^|, while the in- 
flationary paradigm Q has, on the contrary, been widely 
confirmed and is now accepted as part of the "standard" 
model of cosmology. 

Some of the primordial cosmological backgrounds pro- 
posed in the literature are quantum cosmological models 
which share the attractive properties of exhibiting neither 
singularities nor horizons 0, Q, @ , leading the universe 
evolution through a bouncing phase due to quantum ef- 
fects, and a contracting phase before the bounce. They 
constitute an example of a bouncing model, a possibil- 
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ity that has attracted the attention of many authors 
without the presence of any phantom field. These fea- 
tures of the background introduce a new picture for the 
evolution of cosmological perturbations: vacuum initial 
conditions may now be imposed when the Universe was 
very big and almost flat, and effects due to the contract- 
ing and bouncing phases, which are not present in the 
standard background cosmological model, may change 
the subsequent evolution of perturbations in the expand- 
ing phase. Hence, it is quite important to study the evo- 
lution of perturbations in these quantum backgrounds to 
confront them with the data. 

The present paper is the fourth of a series 0, 0, IT^ 
where the theory of cosmological perturbations is ob- 
tained and simplified without assuming any dynamics 
satisfied by the background. This is a necessary pre- 
requisite if one wants to study the propagation of per- 
turbations on a quantized background. The usual the- 
ory of cosmological perturbations with their simple equa- 
tions relies essentially on the assumptions that the 
background is described by pure classical General Rela- 
tivity (GR), while the perturbations thereof stem from 
quantum fluctuations. It is a semiclassical approach, 
where the background is classical and the perturbations 
arc quantized. A full quantum treatment of both back- 
ground and perturbations has already been constructed 
in Ref. |13(, but rather complicated equations were ob- 
tained, even at first order, due to the fact that the back- 
ground does not satisfy classical Einstein's equations. 
In Refs. [ToL IT^ . we have managed to put these com- 
plicated equations in simple forms, very similar to the 
equations obtained in Ref. through the implementa- 
tion of canonical transformations and redefinitions of the 
lapse function without ever using the background classi- 
cal equations. These expressions happen to become iden- 
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tical to those of Q when the background behaves as a 
pre-determined function of time, which is perfectly con- 
sistent with the idea of quantization if we work with an 
ontological interpretation of quantum mechanics |14| . In 
Ref. [ll|], these results have been applied to obtain the 
possible power spectra of tensor perturbations in different 
quantum cosmological models using an hydrodynamical 
description. The aim of this paper is to apply the results 
of Ref. [TjI to obtain the power spectra of scalar pertur- 
bations in such quantum models and confront the results 
with observational data [3j. 

The paper is organized as follows. In Sec. ^ we sum- 
marize the results of Refs. [l'2t obtaining the simple equa- 
tions which govern the dynamics of quantum perturba- 
tions in the quantum backgrounds of Refs. 0, 13- in 
Sec. IIIII we obtain the spectral index for long wavelengths 
of scalar perturbations in these quantum backgrounds, 
and in Sec. II VI we confirm these results numerically, also 
obtaining their amplitude and constraining the free pa- 
rameters of the theory with observational data. We end 
in Sec. V with conclusions and discussions. 



II. QUANTIZATION OF THE BACKGROUND 
AND PERTURBATIONS 

The action we shall begin with is that of General Rela- 
tivity (GRVwith a perfect fluid, the latter being described 
as in Ref. l] 1 , i.e. 



S = S, 



GR i - Sfluid — g^ 2 



—gRd x — 



(1) 

where £ p] = (8nG N /3) 1 / 2 is the Planck length in natural 
units (H = c = 1), e is the perfect fluid energy density 
whose pressure p is provided by the relation p = toe, lo 
being a nonvanishing constant. 

Let the geometry of spacetime be given by 



9^v 



q (0) 



(2) 



where gf?} represents a homogeneous and isotropic cos- 
mological background, 

ds 2 = g^dx^dx" = N 2 (t)dt 2 - a 2 (t)5 ij dx i dx j , (3) 

where we are restricted to a flat spatial metric, and 
the h^u represents linear scalar perturbations around it 
which we decompose into 



h m = 2N 2 <f> 

hm = -NaB^ 

h i3 = 2a 2 (-07ij - E t ij). 



(4) 



1 One can also use the formalism due to Schutz ll "il . obtaining the 
same results. 



Substituting Eqs. (@J and © into the Einstein-Hilbert 
action Q), performing Legendre and canonical transfor- 
mations, redefining N with terms which do not alter the 
equations of motion up to first order, all this without ever 
using the background equations of motion, the Hamilto- 
nian up to second order is simplified to (see Ref. for 
details) 



H = N 



where 



and 



H, 



(o) 



H, 



(2) 



(0) 



+A N P N + j d 3 X(/)TT^+ I d 3 xA07T0, 

(5) 
(6) 



)2p2 



4aV 



3cj 



a: 



(2) 



1 

2~^ 3 



(7) 



The quantities N, (j), Ajy and play the role of Lagrange 



multipliers of the constraints H { 



(0) 



H, 



(2) 



0, 7r^ sa 0, 



Pn ~ 0, and tt^ w 0, respectively. The momenta P a , n^, 
Ti" V" and Pt are conjugate respectively to a, <j>, ip, N, T, 
this last variable playing the role of time. 

The variable v is related with the gauge invariant 
Bardeen potential $ (see Ref. 0) through 



l)eo / v 
ui \ a 



(8) 



which coincides with equation (12.8) of Ref. [| relating 
v and $ 



where 



7 =1 



3 £ P1 Hj 
2 c? 



n 2 ' 



(;)'■ 



(9) 



(10) 



and c s (c 2 = dpo/deo = Po/ e o) ^ s ^ ne velocity at which 
density perturbations propagates, and when the classical 
equations of motion, that can be recast in the form 



f-a+Pa = eo(l +w) 



2H 2 j 
Mia 2 



are used. 

The above quantities have been redefined in order to 
be dimensionless. For instance, the physical scale factor 
flphys can be obtained from the dimensionless a present 
in © through a p h ys = £ P1 a/VV, where V is the comoving 
volume of the background spacelike hypersurfaces, which 
we suppose to be compact. The constraint Hq +Hq is 
the one which generates the dynamics, yielding the cor- 
rect Einstein equations both at zeroth and first order in 
the perturbations, as can be checked explicitly. The oth- 
ers imply that N, <f>, and ip are not relevant. The unique 
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perturbed degree of freedom is V, as obtained in Ref. 0. 
We would like to emphasize again that in order to obtain 
the above results, no assumption has been made about 
the background dynamics: Hamiltonian (JSj is ready to 
be applied in the quantization procedure. 

In the Dirac quantization procedure, the first 
class constraints must annihilate the wave functional 
* [TV, a, <j>(x% »(**), T] , yielding 



9 T 6 T 5 T 

ON 8<j> Sip 



= 0. 



(11) 



The first three equations impose that the wave functional 
does not depend on N, <p and ip: as mentioned above, N 
and 4> are, respectively, the homogeneous and inhomo- 
geneous parts of the total lapse function, which are just 
Lagrange multipliers of constraints, and tp has been sub- 
stituted by v{x l ), the unique degree of freedom of scalar 
perturbations, as expected. 

As Pt appears linearly in H, and making the gauge 



choice N 



one can interpret the T variable as a 



time parameter. Hence, the equation 

= 

assumes the Schrodinger form 



(12) 



• 9 T 

dT 



,(3w-l)/2 



d_ 

da 



,3u-l 



<Tx 



hi 



(3w-l)/2_ 



d xv l v^, 



(13) 



where we have chosen the factor ordering in a in order to 
yield a covariant Schrodinger equation under field redef- 
initions. 

If one makes the ansatz 



*[a,v,T} = $ (0) (a,T)$ (2) [a,v,T], 



(14) 



where ty/ \(a,T) satisfies the equation 



d 



i—* {0) (a,T) 



j(3w-l)/2 Q 



fl (3^-l)/2_^_ 

da 



4 da 

then we obtain for ^ {2) [a, v, T] the equation 



*(o)(a,T), 
(15) 



^*(2)K^ T ) = 2 



/ d3 ^*( 2 )K^ T ) 



,(3^+1) 



d 3 xw' l w^*( 2 )(a, v,T). 

(16) 



Terms involving ^(2) (a, t>, T) can be neglected in 
Eq. I|15l) . either through a judicious choice of the a de- 
pendence of ^(2) |ll|. or because quantum perturba- 
tions initiated in a vacuum quantum state should not 
contribute to it 0. In another perspective, using 
the Bohmian approach [T^j . one can write ^[a,v,T] = 
\E , (o)(a,T)\E , (2)[ w j^ 1 ] an d consider a(T) appearing in the 
equation for ^f^[v,T] as a prescribed function of time, 
the quantum Bohmian trajectory, obtained from the ze- 
roth order equation for $ rm(a,T). 

Going on with the ontological Bohm-de Broglie inter- 
pretation of quantum mechanics, where quantum trajec- 
tories can be defined, Eq. I|16fl can be further simplified 
if one uses Eq. (|15fl to obtain background quantum 
Bohmian trajectories a(T) as in Refs. [UlIlJ- This can 
be done as follows: we change variables to 

obtaining the simple equation 

.a* (0) (o,r) _ ia 2 * (0) (o,T) 



dT 



d\ 2 



(17) 



This is just the time reversed Schrodinger equation for 
a one dimensional free particle constrained to the posi- 
tive axis. As a and x are positive, solutions which have 
unitary evolution must satisfy the condition 



(o)- 



(o) 



dx 



(0)- 



dx 



= 



(18) 



(see Ref. Q for details). We will choose the initial nor- 
malized wave function 



\T ir) y \ T 



(19) 



(init) 
(0) 



where T is an arbitrary constant. The Gaussian \1/ 
satisfies condition (|18|l . 

Using the propagator procedure of Refs. we ob- 

tain the wave solution for all times in terms of a: 



*(0)(a,T) 



8T 


1/4 


exp 





-4T a 



3(1- 



9(t 2 + t 2 )(i-^) 2 



exp 



4Ta 3 ( 1 " 



9(T 2 + T 2 )(l-c) 2 



arctan 





7T 


(?) 


~ 4 



(20) 



Due to the chosen factor ordering, the probability den- 



sity p(a, T) has a non trivial measure and it is given by 
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FIG. 1: Time evolution of the scalar mode function for the equation of state u — 0.1 in the one-fluid model of the bounce. The 
left panel shows the full time evolution which was computed, i.e., the function S(x), while the right panel shows v(x) itself, 
both plots having k — 10~ 3 . For x < 0, there are oscillations only in the real and imaginary parts of the mode, so the amplitude 
shown is a non oscillating function of time. It however acquires an oscillating piece after the bounce has taken place. 



p(a,T)=a^y^ (0) (a,T)\ 
coming from Eq. 1(17(1 reads 



Its continuity equation 



dp 
df 



d_ 

da 



-P 



= 0, 



2 da' 

which implies in the Bohm interpretation that 

ffl (3«-l) Q S 



a = — - 



da ' 



(21) 



(22) 



in accordance with the classical relations d = {a, H} = 
_ a (3c^i)p a / 2 and p a = dS / da _ 

Inserting the phase of 1(20(1 into Eq. 1)22(1 ■ we obtain the 
Bohmian quantum trajectory for the scale factor: 



a{T) = a 



(23) 



Note that this solution has no singularities and tends to 
the classical solution when T — ► ±oo. Remember that we 
are in the gauge N — a 3uJ , and T is related to conformal 
time through 



NdT = adrj 



dr] 



[a(T)] 3w_1 dT. 



(24) 



The solution (|23|) can be obtained for other initial wave 
functions (see Ref. 0)- 

The Bohmian quantum trajectory a(T) can be used 
in Eq. 116(1 . Indeed, since one can view a(T) as a func- 
tion of T, it is possible to implement the time dependent 
canonical transformation generated by 



U = exp 



x exp 



d 3 a 



2a 



mil 

a 



(25) 
(26) 



As a(T) is a given quantum trajectory coming from 
Eq. 1(17(1. Eq. ((25|) must be viewed as the generator of 
a time dependent canonical transformation to Eq. 1(17(1 . 
It yields, in terms of conformal time, the equation for 
®(2)[v,a(Ti),rj\ 



9* 



(2) 



dr} 



d A x 



2Sv 2 



-V iW 



2a 



(2)- 



(27) 

This is the most simple form of the Schrodinger equa- 
tion which governs scalar perturbations in a quantum 
minisuperspace model with fluid matter source. The cor- 
responding time evolution equation for the operator v in 
the Heisenberg picture is given by 



UJV 



0. 



(28) 



where a prime means derivative with respect to confor- 
mal time. In terms of the normal modes Wfc, the above 
equation reads 



uk \Vk = 0. 

a 



(29) 



These equations have the same form as the equations for 
scalar perturbations obtained in Ref. jlj. This is quite 
natural since for a single fluid with constant equation of 
state ui, the pump function z" / z obtained in Ref. [l| is 
exactly equal to a" /a obtained here. The difference is 
that the function a(rf) is no longer a classical solution 
of the background equations but a quantum Bohmian 
trajectory of the quantized background, which may lead 
to different power spectra. 
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III. SPECTRUM OF SCALAR 
PERTURBATIONS IN QUANTUM 
COSMOLOGICAL MODELS 

Having obtained in the previous section the propaga- 
tion equation for the full quantum scalar modes, namely 
Eq. 129fl . in the Bohmian picture with the scale factor 
assuming the form J22J, it is our goal now to solve this 
equation in order to obtain the scalar perturbation power 
spectrum as predicted by such models. In this section we 
obtain the analytical result for long wavelength perturba- 
tions through a matching procedure, while in the follow- 
ing we confirm our current findings by getting numerical 
solutions. 

We shall begin with the asymptotic behaviors. When 
|T| ^S> \Tq\, far from the bounce, one can write Eq. i|29|) 
as 



v + 



2(3uj - 1) 



(30) 



(31) 



(1 + 3w) V 
whose solution is 

with 

3(l-o;) 
2(3w + 1)' 



ci and c 2 being two constants depending on the wave- 
length, iJ(b2) being Hankel functions, and k = ^/ujk. 

This solution applies asymptotically, where one can im- 
pose vacuum initial conditions as in [l| 



e ikr] 



(32) 



377 -<f("+i) 



which implies that 

ci = and c 2 



The solution can also be expanded in powers of k 2 
according to the formal solution Q] 



v 
a 



A 2 (k) 



1 - Luk 2 
' dfj 



a 2 iff) 



a 2 (77) d?y 



dr, 
72" 



* 2,- r # 

a drj 



(33) 



up to order ©(fc-?- 4 ) terms. In Eq. Q33JI . the coefficients 
Ai and A 2 are two constants depending only on the 
wavenumber k through the initial conditions. Although 
this form is particularly valid as long as Lok 2 <C a" /a, 
i.e. when the mode is below its potential, Eq. (|33|l should 
formally apply for all times. In the matching region 



Lok 2 w a" /a, the 0{k 2 ) terms may give contributions to 
the amplitude, but they do not alter the fc-dependencc 
of the power spectrum. 

For the solution (|23[) . the leading order of the solution 
(|3*3^l reads 



Ai + A 2 T a 3 Q UJ ^ ^arctanrr 
Ax-A^oa 3 ^-, 



(34) 



where 



T 

x = —. 

To 

In the last step we have taken the limit x — > — 00, and 
the constant 7r/2 was introduced in order for A\ represent 
the constant mode when it enters the potential. 

Propagating this solution to the other side of the 
bounce, in the expanding epoch, i.e. the limit for x — > 
+00, yields 



v 
a 



Ai 



1 



7r a, 



3(w-l) 



T A 2 



= [A 1 +Tral^- 1) T A 2 



3(w-l) 

V 



T A 2 . (35) 



Note that there is a mixing in the constant part of the 
mode when it passes through the bounce. Hence, in these 
types of bouncing models, the bounce has important con- 
sequences for the final power spectrum. 

In order to find the fc-dependence of A\ and A 2l we 
match v and v' in Eqs. H34|) and (|31|l . obtaining, to lead- 
ing order, 



A l = 



A 2 = 



3cj- 1 
3a (uj - 1) 



C 



T (3D 



_ 3(l-m) 
\ 2(3to+l) 



2a, 



l+3w 




30(l-u)T o 



C-aD 



_ 3(^-1) 
£ \ 2(3u + l) 



from which stem the spectral behaviors 



3(l-i^) 



3(^-1) 

Ar, cx fc 2 < 3 -+ 1 > . 



(36) 
(37) 

(38) 



with 



9(1- w) 2 
2 II — 3^1 



-l/(l+3w) 



= 



9(1-0;)' 
211 - 3w| 



C = VTo~a^- 1)/2 c 2 ^k^Hl 2 Hk^ 



(39) 



and 



n _ (i-3w)/2„-i/2C2 J Hi >(kr) u ) 

T) n — 

1 { V k Vm 



+ Jkr ]M Hi 2 \(k VM ) - Hi%(kr, u ) ,(40) 



6 



with 



and 



V2|l-3o;| 
l + 3w 



T c% 



3u)-l 



(41) 



The coefficients A\ and each contain a power-law 
behavior in k. Because < u> < 1, the power in A 2 
[Eq. is negative definite and that in A\ [Eq. (J2HJ)] 

is positive definite. Therefore, A2 is the dominant mode, 
while A\ provides the sub-dominant mode. 

The relation between the Bardeen potential $ and v is 
given by Eq. JHJ). Using Eqs. C£9 and J23J in order 
to change variables to x, Eq. JSJ leads to 



(l+3u»)/2 



<I>a 



A 2 (k) 



a 
-uik 







1 — ujk 2 


f dxa 2 J 







^ J X T a 2 (x)dx 
+ 0(k 4 ) 



T 2 



a 2 (l + x) 



:dx 



2 MQ 
a 2 



Too 2 (x) dr] 



A 2 (k) 



l-Wk 2 ^ 

a, 



'2 r-x 

y / dxa 2 (arctanx + 7r 
J 



/2) 



(42) 



where the constant 7r/2 was introduced when performing 
the last integral in order for A\ to represent the constant 
mode at large negative values of x when is entering 
the potential, as before. At large positive values of x, 
when Vk is leaving the potential and x oc 7y 3 ( 1_w )/( 1+3w ), 
the constant mode of like v, mixes A\ with Ai. In 
this region, taking into account that Ai dominates over 
Ai, we obtain: 



const. 



1 



(5+3^)7(1+3^) 



The power spectrum 



■K 2 



then reads 



and we get 



oc k n s- 



n a = 1 



I2u 
1 + 3lj' 



(43) 



(44) 



(45) 



(46) 



For gravitational waves (see Ref. [llj for details), the 
equation for the modes [i = h/a reads 



2/C ]a = 0, 

a 



(47) 



yielding the power spectrum for long wavelengths: 

2k 3 11, 2 



■K 2 



In Ref. [Tl|, we have obtained 
V h oc ¥• 



with, as for the scalar modes, 



12cj 



(48) 



(49) 



(50) 



Note that in the limit ui — > (dust) we obtain a scale 
invariant spectrum for both tensor and scalar perturba- 
tions. This result will be confirmed in the next section 
through numerical calculations which will also give the 
amplitudes. However, it is not necessary that the fluid 
that dominates during the bounce be dust. The depen- 
dencies on k of A\ and A2 are obtained far from the 
bounce, and they should not change in a transition, say, 
from matter to radiation domination in the contraction 
phase, or during the bounce. The effect of the bounce 
is essentially to mix these two coefficients in order for 
the constant mode to acquire the scale invariant piece. 
Hence, the bounce may be dominated by another fluid, 
like radiation. If while entering the potential the fluid 
that dominates is dust-like, then the spectrum should be 
almost scale invariant. Note also that since we assume 
an ordinary matter fluid, the equation of state is positive 
definite, being, in the most pressure-less case, obtained 
as a result of the quantum non vanishing mean square 
velocity. Thus, we do expect a blue spectrum. 



IV. NUMERICAL RESULTS 

In this section we will confirm the spectral indices 
of scalar perturbations obtained above, and obtain the 
amplitude of the scale invariant mode. The dynamical 
mode equation is expressible in terms of the function 
S = a^^ 1 ^^) v / \/Tq (the constant VTo being introduced 
in order for S to be dimensionless), namely [TlJ 



S- 



P(l 



1 



(1 



5 = 0, (51) 



with, in this latter case, x = T/Tq and k = c s /cTo/oJ u 
We can apply the vacuum initial conditions 



^ic s krj 



(52) 



with the sound velocity c s = ^fu> a constant. It is clear 
here that one must insist upon not having u> — in order 
to be able to put these initial conditions. Again, anyway, 
the sound velocity, even in a matter dominated phase, is 
not expected to be vanishing identically. 
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FIG. 2: Rescaled power spectra for scalar and tensor pertur- 
bations as functions of time for uo — 0.01 and two different 
values of k. It is clear from the figure that not only both 
spectra reach a constant mode, but also that this mode does 
behave as indicated in Eqs. 14611 and 15011 . It is purely inci- 
dental that the actual constant value of both modes are very 
close for that particular value of u. The constant values ob- 
tained in this figure are the one used to derive the spectrum 
below. In this figure and the following, the value of n s used 
to rescale $ is the one derived in Eq. 1461 . thus proving the 
validity of the analytic calculation. 



FIG. 3: Rescaled power spectra for u = 8x 10 -4 , correspond- 
ing to the conservative maximum bound on the deviation from 
a scale invariant spectrum n s = 1.01, as function of k. The 
scalar spectrum V<s> is the full line, while the dashed line is 
the gravitational wave spectrum Vh- Also shown is the ratio 
T/S (dotted); in this case, the T/S ~ 5.2 x 10~ 3 , i.e. almost 
two orders of magnitude below the current limit. This case 
has a typical bounce length-scale of Lo ~ 1.47 x 10 3 ^ pl . The 
amplitude of the modes is obtained as the constant part of 
Fig.H 



Defining 



L,Q — TqQq^ 



(53) 



the curvature scale at the bounce (the characteristic 
bounce length-scale), namely, Lq oc 1/^/Rq where Rq is 
the scalar curvature at the bounce, from which one can 
write 



«o 



-L 



\ bounce ' 
phys 



(54) 



one obtains the scalar perturbation density spectrum as 
a function of time through 



2( " +1) 7 mf n 



7T 2 V^(1-W)) 2 fc 

where 

f(x) = (l + a 
while the tensor power spectrum is 
2fc 3 \v\ 2 ft 

V h = 



(55) 



S, 



n 2 1 + x 2 \L 
in which the rescaled function v, defined through [Tll | 

v — ■ , 



(56) 



satisfies the same dynamical equation 1)510 with c s 
with fx subject to initial condition 



^—ikr) 



l, 



(57) 



From the above defined spectra, one reads the ampli- 
tudes 



and 



A 



2 



lilt) 



(58) 



(59) 



where we assume the classical relation between $ and the 
curvature perturbation £ through 



c 



5 + 3m 
3(1+ «) 



3lu 



3(l + a;) 



V<£, 



(60) 



to obtain the observed spectrum. Since both spectra are 
identical power laws, and indeed almost scale invariant 
power laws, the tensor-to-scalar (T/S) ratio, defined by 
the CMB multipoles C t at i = 10 as 



(T) 



T ^ C w 



S 



c 



(S) 
10 



A 2 



(61) 



can easily be computed (see, e.g., and references 
therein). In Eq. I|61|) . the function T depends entirely 
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on the background quantities such as the total energy 
density relative to the critical density, i.e. f2, among oth- 
ers. It does not depend on primordial physics parameters, 
supposed to be included in the amplitudes of the spectra; 
in other words, it propagates the predicted primordial 
spectra through the different epochs (radiation, matter 
and cosmological constant dominated) whose characteris- 
tics are fixed by observational data. For Eq. (|61|) to hold 
true, both scalar and tensor spectra ought to be power 
laws, as they are in our model. 

Eq. (|61l) permits to compare primordial cosmological 
effects such as the ones derived here with current ob- 
servations. This means that our model must somehow 
be connected with the observed universe, which we take 
to be the so-called concordance model (the one having 
a cosmological constant accounting for fl\ ~ 0.7 of the 
total density). The relevant value of the function T can 
be obtained from the slow-roll result: in this case, on 
finds [l6j A^/A^ ~ 16e, with e the usual slow-roll pa- 
rameter, whereas the consistency equation demands that 
C^o /CIq* ~ lOe. For the concordance model, one thus 
has T ~ 10/16 ~ 0.62. This is the value we shall use to 
compare our models with observations. 

We have solved these equations numerically and ob- 
tained the values of the free parameters which best fit 
the data. First, we assume a spectral index limited 
by n s < 1.01, an admittedly conservative constraint, 
which, given Eq. |@5J|, provide the already severe bound 
on the equation of state ui < 8 x 10~ 4 . Constraining 
= 2.08 x 10~ 10 then implies the characteristic bounce 
length-scale Lq to be 

L > 15004,, 

a value consistent with our use of quantum cosmology: 
it is indeed in the kind of distance scale ranges that one 
expect quantum effects to be of some relevance, while at 
the same time the Wheeler-de Witt equation to be valid 
without being possibly spoiled by some discrete nature 
of geometry coming from loop quantum gravity, and/or 
stringy effects. The tensor-to-scalar ratio is, in this case, 
much lower than the current WMAP3 constraint 01 
T/S < 0.21 as we obtained T/S ~ 5.2 x 10~ 3 . It is 
interesting to note that when the T/S constraint is al- 
most reached, this gives a numerical estimate for lu which 
is u> ~ 8.5 x 10~ 2 , and hence a spectral index n s ~ 1.81, 
i.e. much larger than the current constraint. In this case 
also, one finds Lq ~ 350^ pl . Hence, it is not necessary to 
explore these possibilities. 

Let us briefly discuss the validity of the model dur- 
ing the bounce epoch. It is indeed not enough that the 
Wheeler-de Witt equation applies in the regime one con- 
siders, but we should also make clear that the backre- 
action is not going to dominate at the bounce. That 
would for instance be the case if the Bardeen potential 
grows large enough. In other words, we must ensure that 
\$(x)\ <C 1 at all times, a condition which, in view of 
Fig. 121 might not be easily satisfied. Indeed, in a way al- 
most independent of u>, one finds that the spectrum V®, 



once normalized to the CMB for late times, is of order 
unity (depending only very little on the values of k and 
lu) at its maximum, close to the bounce. This is not so 
much of a problem since what needs to be small is the 
Bardeen potential itself, and its Fourier mode reads 

- - i^f) ■ m 

having an explicit dependence in ao- As this latter pa- 
rameter is not constrained (in fact, it must be such that 
large physical scales now must be much greater than L 
at the bounce in order for our approximation k -C 1 be 
valid all along, see Eq. i|54|)L one can assume it is suffi- 
ciently large to yield <C 1 at all time and for all k of 
cosmological relevance. Note that ao could in principle 
be fixed by the normalization of the scale factor at the 
bounce in terms of, say, the Hubble constant today, to 
fix the scales A = a/k to be in the observable range now. 
This will be done when we examine in a future publica- 
tion a more elaborated model containing not only dust 
but also radiation, as discussed in the end of Sec. Ill, 
in such a way that it can be connected to a radiation 
dominated expansion phase before nucleosynthesis. 

V. CONCLUSIONS 

Using quantum cosmology and the Bohm interpreta- 
tion, thanks to which one can define trajectories and a 
scale factor evolution with time, we have obtained a sim- 
ple model whose scalar and tensor perturbations can be 
made arbitrary close to scale invariance. The model con- 
sists of a classically contracting single dust perfect fluid 
in which the big crunch is avoided through quantum ef- 
fects in the geometry described by the Wheeler-DeWitt 
equation, turning the Universe evolution to an expanding 
phase, which soon becomes classical again. This transi- 
tion is smoothly described by a Bohmian quantum tra- 
jectory containing a bounce. Perturbations begin in a 
vacuum state during the contraction epoch, when the 
Universe was very large and almost flat, and are subse- 
quently evolved in a fully quantum way through its all 
history. Hence, we have presented a nonsingular model 
without horizons where perturbations of quantum me- 
chanical origin can be described all along and generate 
structures in the Universe. 

In order to match the CMB normalization, we find that 
the characteristic length at the bounce must be of the or- 
der of a few thousands Planck lengths, thereby making 
the model fully consistent. Indeed, we should like to em- 
phasize that one would expect precisely the Wheeler-de 
Witt equation to be valid and important in this regime (a 
sufficiently low scale for quantum gravity effects be im- 
portant, but not so low in order for not being affected by 
string and/or loop quantum gravity effects). Moreover, it 
predicts a slightly blue spectrum, which may be seen as a 
drawback of the model, although this point still deserves 
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further experimental clarification, say by the Planck mis- 
sion. However, in more realistic and elaborated models, 
other fluids must be considered, like radiation and dark 
energy. It seems that adding these fluids will not spoil 
our results as long as only a single fluid dominates at 
the bounce, and that dust dominates when the scales of 
cosmological interest become greater than the curvature 
scale in the classical contracting phase. If in that period 
dark energy has some effect imposing a slightly negative 
effective w, then one could even obtain a slightly red 
spectrum instead of the blue one derived here. Besides, 
ironically coming back to one of the original motivation 
for inflation [l^ , we remark on Fig. [3] that the gravita- 
tional potential shows a net increase for large values of 
k after the bounce has taken place. That could initiate 
non-linear growth of the primordial perturbations at very 
small scales (provided the values of To and oo are chosen 
conveniently). In turn, these would lead to the formation 
of primordial black hole whose decay could initiate a ra- 
diation dominated phase, as needed. These elaborations 
should also consider other classical cosmological puzzles, 



namely the flatness and remnants (e.g. monopoles) prob- 
lems, and baryogenesis. These issues will be considered 
in future works. 

Finally, we have obtained that once the constraint on 
the scalar index is taken into account, the tensor-to-scalar 
ratio which follows is predicted to be small. Thus, this 
category of model is falsifiable. 
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